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Abstract 

The complete 1-loop computation of the processes (xfXj 77 > 99) > f° r an y P a i r °f the 
four MSSM neutralinos, has been performed for an arbitrary cm energy. As a first ap- 
plication suitable for Dark Matter (DM) searches, the neutralino-neutralino annihilation 
is studied at relative velocities v/c ~ 10~ 3 describing the present DM distribution in the 
galactic halo, and v/c~ 0.5 determining neutralino relic density contributions. Numerical 
results are presented for 31 MSSM benchmark models indicating considerable sensitivity 
to the input parameters. Numerical codes are released, which may be used for the com- 
putation of the annihilation of any two neutralinos to two photons or two gluons at the 
aforementioned v/c values. In the near future, we intend to complement them with codes 
for the inverse process 77 — > XiXj> observable at the future high energy Linear Colliders. 

PACS numbers: 12.15.-y, 14.80.Ly, 95.35+d 
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The nature of the Dark Matter (DM), which is here assumed to be possibly sensitive 
only to the usual gravitational and weak forces PjJ |2], is still open. At present, a most 
obvious candidate for such matter in the context of R-parity conserving Supersymmet- 
ric (SUSY) theories, should be the Lightest Supersymmetric Particle (LSP) 3J. In the 
minimal Supersymmetric models (MSSM), this is most often identified with the lightest 
neutralino Xi |UE]- Within MSSM there exist cases though, where the LSP is the purely 
gravitationally interacting gravitino, while the lightest neutralino may then appear as the 
next to lightest supersymmetric particle (NLSP) jB]. 

In any case, depending on the MSSM parameters, one or more neutralinos may sizably 
contribute to DM [TJI2]- Such DM could be investigated though direct or indirect methods. 
Direct detection could occur through e.g. the observation of the nucleus recoil induced 
by the x^-nucleus scattering [7J. Experimental setups are looking for such events jH], and 
there even exist claims that a signal may have already been observed jHj- 

Indirect detection refers to searches of anomalies in the spectrum of the photons, 
neutrinos, positrons or antiprotons, emitted from cosmologically nearby sites, where the 
neutralino concentration might be enhanced PjJ |2] . Both types of searches are however 
difficult, because of uncertainties in making precise predictions for the dark matter density 
distribution in these sites ^U], and the strong sensitivity of the observable quantities on 
the SUSY parameters. 

However, as noticed several years ago [TTJ H2] , the process 1 XiXi ~~ > 77 ( as we U as 
X1X1 ~ * should be particularly interesting for Dark Matter detection, since the 

neutralino pair annihilation is essentially taking place at rest, producing sharp monochro- 
matic photons. Provided therefore that the cross section cr(x?Xi — > 77) an d the neutralino 
concentration in the galactic halo are not very small, the signal to be searched for would 
be sharp photons coming from the central region of the Galaxy, with energies in the 
hundred-GeV region. 

The companion process X1X1 99 ma J a lso be important in some regions of the 
parameter space, where its sufficiently frequent occurrence may modify the analyses of 
the neutrino and antiproton signals ^3]. Furthermore, if the second neutralino x\ turns 
out to be very close to £?, processes like (xiX.2 -> 77, 99) and (xlxl -> 77, 99) may 
also be interesting; particularly the gluonic mode may then play a role in determining the 
relic neutralino density and the rates of neutrino and antiproton production [T3] . 

These statements sufficiently motivate our interest in the general processes x^X^ 
77, 99 1 {hi = 1 — 4 ). The lowest order non- vanishing contribution to them arise from one 
loop diagrams involving standard and supersymmetric particle exchanges; i.e. sleptons, 
squarks, gauginos and additional Higgses, which sensitively depend on the precise SUSY 
particle spectrum, couplings and mixings. Thus, DM considerations of the processes 
(XiXj ~^ 77, 99)i an d studies of the inverse processes 77 — > XiX.'j m a Linear Collider 
(LC), and gg —>■ x^x'j a t LHC, may provide very strong constraints on the various MSSM 
models. 

Previous computations have only addressed the diagonal annihilation XiX^ ^77, 99 
1 Xi denotes the lowest neutralino. 
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of two identical neutralinos at vanishing relative velocity [TTjrE^ rEJ]. Such approximations 
may of course be adequate for neutralinos annihilating at present in the galactic halo (or 
the Sun or Earth), where their average relative velocity is expected to be v/c ~ 10~ 3 
[I]. But this low energy limit may not be generally adequate for neutralino relic density 
computations, which is determined by neutralino annihilation at a time when their average 
relative velocity prevailing in Cosmos was v/c~ 0.5 fTl I14j. 

We have, therefore, made a complete one loop computation of 11 ,99, m a 

general MSSM model with real parameters, expressing all helicity amplitudes in terms 
of the Passarino-Veltman functions, for general incoming and outgoing momenta. The 
derived expressions, which should be useful in both dark matter and LC and LHC studies, 
are quite complicated though, so that the only efficient way to present them is through 
numerical codes. Details for both types of processes and applications to collider physics 
searches, will be given in a separate extended paper [To] . 

In the present paper we restrict to studies of the xfXj — y 11 and XiXj "~ > 99 an- 
nihilations in the framework of DM searches, where the average relative velocity of 
the two neutralinos is small. Two numerical codes, PLATONdml and PLATONdmg 
are released which compute the neutralino annihilation cross sections times v^; i.e. 
(vij/c) ■ a(x iX j -> 77) and (v y /c) • cr(x°X° -> 99) respectively, at v^/c ~ 10" 3 , as it 
would be expected for galactic halo neutralinos at present [TTJ O E] . In addition, a third 
code PLATONgrel, is also released, which calculates (%/c) ■o-(XiX < j ~~ > 99) at v^/c ~ 0.5, 
which should be useful for neutralino relic density computations. The outputs of all codes 
are in fb. 



Procedure . 

The generic structure of the one-loop diagrams for Xi Xj — > 11 m the 't-Hooft-Feynman 
gauge, is depicted in FigJTJ where the types of particles running clockwise inside each 
loop are: box (a), (fSSS), (fWWW), (fSSW), (fSWW), (fSWS), (fWWS); box 
(b), (Sfff), (Wfff); box (c) (SffS), (WffW), (SffW), (WffS); initial triangle (d), 
(SfS), (WfW); final triangle (e), (WWW), (SSS), SWW , SSW; final triangle (f), 
(///); and bubbles (g) (SS), (WW); (h). WS. By S we denote the scalar exchanges 
(Higgs, Goldstone and sfermions); and by / the fermionic ones (leptons, quarks and inos). 
Bubbles (g) and final triangles (e) and (f), are connected to the initial XiXj state through 
an intermediate Z, or neutral Higgs or Goldstone boson h°, H°, A , G°. For the x°Xj — * 99 
process, only the subset of the diagrams in FigJU involving quark and squark loops, is 
needed. We also note that the XiXj-antisymmetry, and the Bose symmetries of the 77 
and gg states, provide non-trivial tests for the helicity amplitudes of XiX^j ~^ 11, 99, 
which we have checked to be satisfied [15]. 

General features for small relative velocities. 

When the relative velocity among the neutralino pairs is small, their cm energy to 
order in (t^/c) 4 , is 
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For very small values of 2 Uy, the angular distribution of da{x^,X°j —> 77, 99) /d cos 9 is 
flat, and %cr(x°x° — > 77, gg) is a smooth function of u?-. 

Since cx(x°x° ~~ ► 77> fl'fl') are induced by 1-loop processes, the expected order of mag- 
nitude for very small %/c is 



a 4 c 3 x 10 32 cm 3 sec 1 



%*(XiX;->77) - M 2^- M^/TeV) 

,„n~n x a 2 a 2 c 6 x 10~ 30 cm 3 sec~ 1 

ViAxtx^gg) - 7^7" m^m) ' (2) 

where M e // is an effective mass-scale, depending on the masses and mixings of the external 
neutralinos and the particles running inside the loops. For such Vij/c, the two neutralinos 
are predominantly in a 1 So state, so that only helicity amplitudes satisfying Ai = A2 and 
\L\ = fi 2 are allowed. 

The dominant diagrams in this case are the boxes (Sfff), (SffS) (fWWW), (Wfff), 
(WffW), (fWGW), (GffW), (WffG). Triangle and bubble contributions are either 
vanishing or numerically negligible, except for the case where the s-channel exchanged 
neutral Higgs is very close to the sum of the two neutralino masses; compare Fig^-h. In 
such a case, depending on whether x° and xl have the same or opposite CP eigenvalues, 
the main contribution arises from the exchange of either A or 3 if , respectively. 

The magnitude of v^a in (J2J), may thus be enhanced, if relatively light stops, sbottoms, 
and charginos or staus (for the 77 case only) exist, and one of them happens to be almost 
degenerate to one of the incoming neutralinos; or if a neutral Higgs, with an appropriate 
CP eigenvalue, is very close to the cm. neutralino energy, (compare the diagram of 
FigH^-h). But the most important role in determining the v jjO"-magnitude seems to be 
played by the nature of the neutralino contents; mixture percentages of Wino, Higgsinos 
or Bino. The Wino-type neutralinos usually give the largest cross sections, followed up 
by the Higgsinos, while Binos are most often supplying very small cross sections. This 
ordering comes from the dominant role played by the box diagrams involving W bosons 
and the neutralino-chargino-W couplings. For example, it has been emphasized that if the 
lightest neutralino Xi is Higgsino-like and close to a chargino, then M e ff ~ rnw, leading 
to va(xiXi ~^ 77) ~ 10~ 28 cm 3 sec _1 , independently of the neutralino mass; compare (J2J), 
[TT| IT2] . We confirm this result. In fact, we have numerically checked that our exact 
1-loop computation reproduces the results of [TT1 H2] , obtained through a non-relativistic 
treatment of the annihilation cross section of two identical neutralinos at threshold. It 
should also be remarked that (0) is consistent with unitarity in our calculations, provided 
that the neutralino masses are in the TeV-range or below, [TZj. Thus, higher order effects 
can be safely ignored. 

Illustrations for various MSSM models. 

Using the numerical codes PLATONdml, PLATONdmg and PLATONgrel, we obtain the 



2 For all MSSM-models we are aware of, the differential cross section is flat for Wy/c < 0.5. 
3 The lightest h° is below threshold in all models considered. 
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results for Vija(XiXl — > 77) an d Vij&(x'iX ( j ~^ 99) presented in Tables 1, 2 and 3, for an 
extensive set of benchmark SUSY models involving real parameters only. These models 
are used here just in order to indicate the range of possible results we could obtain. 
They all satisfy /i > and have been suggested as being roughly consistent with present 
theoretical and experimental constraints [TH] EH 1201 I2U 122] ■ In constructing these codes, 
the exact 1-loop expressions for XiX^ ~^ 77; 99i have been used. 

In the same Tables, the grand scale parameters of the various benchmarks MSSM 
models are also shown. Conveniently, these models are divided into three classes: 

• Universal m-SUGRA models. Apart from the /x-sign, there exist four additional 
GUT scale parameters in these models, namely (mi/ 2 , m , A , tan/3). They are 
given in the left part of Table 1. The SPiS-models in this Table come from the 
Snowmass set [TH], AD(fgh)x t % from Fig.5 of P3I], and Rosz 2 from Fig. 2 of [20j. 

• Non-universal m-SUGRA models. The grand scale parameters of theses models are 
presented in the upper part of Table 2. In the models AD(fg9) (coming from 
Fig. 9 of [TO]) and SPS6 is allowed to deviate at the grand scale from the 
scalar mass unification condition; while in the CDG- models, the gaugino large scale 
unification is violated |21) . 

• GMSB and AMSB models. The grand scale parameters for the gauge mediated sym- 
metry breaking (GMSB) models (SPS7, SPS8), and the anomaly mediated sym- 
metry breaking (AMSB) models SPS9, are given in Table 3, [T8] . 

Using a public code like e.g. SuSpect j2H], the electroweak scale parameters of the 
above benchmark models may be calculated from their grand scale values; while the also 
needed total widths of the supersymmetric neutral Higgs particles 4 H° and A , are calcu- 
lated using 5 HDECAY [24 . These in turn are used in the in-files of the three PLATON 
codes. 

We next turn to the results for %cr(x°x° — > 77, gg) for the 6 (XiXi), (X1X2) an d 
(X2X2) cases , which appear in Tables 1,2,3. Of course, the (x®, xl — ► 77) case is the most 
interesting one for Dark Matter detection in the halo, but the other modes (especially 
the gluonic ones if they turn out to be appreciable), may play an important role for the 
determination of the value of the DM relic density and the rates of neutrino and antiproton 
production [IB]. Below we discuss these results, in conjunction with the electroweak scale 
masses and couplings of the various models. 

The X1X1 — ¥ 77 mode. 

In all universal mSUGRA models of Table 1, Xi is predominantly a Bino, which implies 
rather small values for v a, in agreement with fU^]. The smallest v a values, in the range 

4 Compare the diagrams in Fig^-h. 

5 In some non-universal SUSY models we are using, an HDECAY version is needed allowing indepen- 
dent values for the gaugino parameters Mi and Mi at the electroweak scale. We thank Abdelhak Djouadi 
for providing us with this code. 

6 Here (xv X2) are the lightest and the next to lightest neutralinos. 
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of 10~ 33 cm 3 sec _1 , appear for the SPS2 models characterized by very heavy sfermions 
and charginos. Such values increase to the 10 -31 cm 3 sec' 1 level, for the R0SZ2, SPS4 
and AD(fg5)2 models, where sfermions and charginos have masses of a few hundreds of 
GeV; and they are further enhanced to a few times 10~ 30 cm 3 sec~ 1 , in models like SPS3i 
or SPS32 where the lightest Xi is close to f\. These enhancements remain nevertheless 
moderate, because the intermediate (xiXi — * t~t + ) -contribution to the X1X1 ~* 77 
amplitudes is depressed by the small r mass. 

The non-universal mSUGRA set of models (Table 2) includes cases like AD(fg9), in 
which the Xi is predominantly of Higgsino-type, leading to va ~ 3 x 10~ 29 cm 3 sec~ 1 
[TT) IT2"]. Much larger values like va ~ 10~ 27 cm 3 sec -1 are reached in the CDG200 and 
CDG011 cases though, in which x\ is predominantly a Wino combined with an appreciable 
Higgsino component. The fact that xt is rather close to Xi m the AD(fg9), CDG200 and 
CDG011 models, also plays a role in enhancing vaXxXi ~^ 77- 

In all other cases of Table 2, Xi is mainly a Bino and va is small. Among them, worth 
mentioning may the the CDG75 model, which displays the very rare property to have a 
CP-odd Xi- 

When it happens that m^ ~ 2m^o, then it is very important to correctly include the 
A width, which we have done by using the code HDECAY [24J. An example of a case 
where the A width is necessary, is CDG^oo- 

In Table 3, we present the results for the GMSB models (SPS7, SPS8), and the 
AMSB models SPS9 [TH]. In (SPS7, SPSS), xi is mainly a Bino and va rather small 
compared to the result in the SPS9 model where Xi is mainly a Wino. In (SPS7, SPS8) 
there exists a near degeneracy of Xi with fi, which enhances the intrinsically small 
Bino result up to the 10 -29 cm 3 sec~ 1 level. In SPS9, which is characterized by a near 
degeneracy of the Wino-like Xi with xt 1 the enhancement is much bigger, leading to 
va ~ 1CT 27 cm 3 sec~ 1 . 

The X1X1 -> 99 mode. 

In Tables 1,2,3 we give results both, for v/c ~ 10~ 3 and v/c~ 0.5. As in the 77 case, the 
largest va values appear for a Wino-like Xii smaller ones for Higgsinos, and very small 
values for the Bino case. On top of this, further enhancements may occur when x\ is 
very close to t\\ or when m^o ~ 2m^ 1 where a considerable sensitivity on v/c may also 
be induced, (see e.g. CDG200 in Table 2). 

As seen from Tables 1,2,3, the va rates vary from 2 x 10~ 31 to 7 x 10~ 29 cm 3 sec _1 
for the universal m-SUGRA models; from 10~ 30 to 10~ 28 cm 3 sec~ 1 for the non-universal 
m-SUGRA models, with the exceptional case of CDG200 giving ~ 10 _26 cm 3 sec _1 ; and 
from 3 x 10~ 31 to 2 x 10 -29 cm 3 sec~ 1 as we move from the AMSB towards the GMSB 
models of Table 3. 

The X1X2 an d X2X2 annihilations. 

In the 77 case, va^XiXV) is very often considerably smaller than va(xiXi), while va^x^X^) 
considerably larger and often reaching the 10~ 27 cm 3 sec~ 1 level. This is usually driven 
by the fact that x\ is mainly a Wino in all models where Xi dominantly a Bino. The 
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occasional degeneracy of X2 with xf, further enhances this value. The only models where 
the X2X2 ra ^ e i s smaller than the X1X1 one > are CDG200, CDGon and the AMSB ones, 
in which the Wino component of Xi is usually considerably larger than that of 

In the gg case, the X1X2 cross sections can be either larger or smaller than the X1X1 
ones, because of the varying situations for the squark and chargino contributions, and 
orthogonalities that tend to appear in certain couplings. In contrast to this, X2X2 is uiost 
often larger than XiXv but the relative magnitudes are often not so large as in the 77 
case. Rates like va{x%xV) ~ 10 -28 cm z sec~ l are often realized. Exceptional cases like 
SPS3i and SPS6 3 give 10~ 26 cm? sec' 1 at v/c ~ 10~ 3 , while SPS6 2 gives a similar value 
at v/c ~ 0.5 at least partly due to A resonance effects. 

Final comments . 

In conclusion, we confirm that fjjCr(x°Xj ~~ > 77; 99) annihilation at Vy/c ~ 10~ 3 , can 
assume a wide range of values. This arises because these processes start at one loop, 
thereby being sensitive to many aspects of the sfermion, ino and Higgs boson masses, and 
the peculiar mass-coincidences appearing in many MSSM cases. 

The results should be useful for Dark Matter searches through neutralino annihilation 
to sharp photons. Roughly, the gamma-ray flux expected in such a case from a direction 
at an angle if) with respect to the direction of the galactic center, is given by [TUj 



dcosif) 



(io-'wv) gfeg 7:7 m , (.->,) 

10 za cm d s 1 



where J (if)) depends on the dark matter density distribution along the line of sight. Its 
magnitude is expected to be of order one at high if); but it increases as if) decreases, and 
in some models it may even reach values up to 10 5 when looking towards the galactic 
center [TU]. On the basis of this, it seems that for va(xiXi ~~ > 77) ^ 10 _29 cm 3 s _1 , the 
observation of a sharp photon line with an energy E 1 precisely corresponding to the mass 
of the LSP Xii should be possible The actual observation of such energetic photons 
would supply an interesting confirmation of Supersymmetry, provided of course that the 
actual observations would allow for acceptable values of the relic density of the lightest 
neutralino. In this respect, the process x^X°j ~ > 7^ may also be interesting, and we are 
planning to extend our study to it. 

The neutralino relic density depends on the total annihilation cross section (summing 
all final states like //, WW, ZZ, Concerning this, we have found that the gluonic 

(sometimes even the photonic channel) may be important in some models, and that occa- 
sionally considerable variations may appear in — > gg), as %/c varies between 
10~ 3 and 0.5. In more detail, the neutralino contribution to cold dark matter (CDM) is 



6 • 10 27 cm 3 s 3 

('-^ann^'rel) 



L CDM 



(4) 



where (<7 ann f re i) is the average total neutralino annihilation cross section multiplied by 
the relative velocity at v rc \/c ~ 0.5. For a model to be consistent with cosmological 
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constraints, we expect that (cr ann t> rel ) ~ 10 _26 cm 3 s _1 . Thus, we note e.g. from Table 2 
that in the CDG 2 oo model, both va(x%Xi ~~ > 99) an d vcr (xiXi ~* 77) are of the order 
10 _27 cm 3 s _1 . A similar situation appears for va(xiXi — > 77) in the AMSB Snowmass 
models [TH] in Table 3. We are tempted to conclude therefore that it would be safer, if 
both these processes are included in the relic density computations. 

Because of the sensitivity of the results to the SUSY models, and in order to contribute 
in the successive stages of the required analysis, we have started the construction of 
a series of numerical codes (PLATON) which will be made public so that anyone can 
run his preferred set of MSSM parameters or explore the parameter space. At present 
the codes PLATONdml allow to compute the annihilation rates (%/c)cr(x°X? — * 77) a ^ 
v/c ~ 10~ 3 , while PLATONdmg and PLATONgrel determine (vij/c)a(xiXj ~^ 99) at 
v/c~ 1CT 3 and v/c ~ 0.5 respectively, always in fb [TH] . 

If the emerging overall picture turns out to be cosmologically consistent, it will give a 
stringent test on SUSY models, particularly if X1X1 ~ > 7^ is a l so observed with the correct 
properties. Particle Physics experiments should then provide further tests by studying 
neutralino production at high energy colliders, through the inverse processes 77 — > xlx'j 
at LC and gg — > XiX^ a ^ LHC. A theoretical study of these processes using the same 
benchmark models as in the present paper is in preparation [To] . 
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7(ei.Pi) xjft) 



(f) 



7( e 2,P2) 




if CO 



7( £ i>Pi) Xi('i) 



7( £ i,Pi) 



(g) 



00 



x,°(/i) 




7(fi>Pi) 



Figure 1: Feynman diagrams for XiXj ~^ 77- Full internal lines denote fermionic ex- 
changes; while broken internal lines denote either scalar or gauge exchanges, except in 
the diagram (h), where the W and Goldstone exchanges are indicated explicitly. 
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Table 1: Universal mSUGRA Models: Input parameters are given at the grand scale, 
and %c(x°x° — > 77, gg) predictions are at %/c ~ 10~ 3 (%/c ~ 0.5). 



Parameters (Dimensions in GeV) 


V 


iAxW - 


-> 77, gg) in 10 29 cm 3 sec 1 








> U 






77 


99 


Model 


m 


mi/2 


A) 


tan (5 


-u -0 
X1X1 


X1X2 


~u ~o 
X2X2 


~o ~o 
X1X1 


-u -0 
X1X2 


-u -0 
X2X2 


CI T~i CI f 

SPS5 


150 


300 


-1000 


5 


0.234 


0.199 


152. 


0.0906 


0.689 


8.40 


















(0.0891) 


(1.52) 


(9.23) 


bFblai 


100 


or n 

250 


-100 


10 


0.443 


0.643 


144. 


0.460 


0.517 


4.40 


















(0.474) 


(0.580) 


(8.00) 


C1 J-} C1 -1 

bPbla2 


140 


350 


-140 


10 


0.291 


0.164 


139. 


0.212 


0.205 


10.6 


















(0.214) 


(0.225) 


(28.7) 


CI J ~> (~M 

6P61a3 


200 


500 


-200 


10 


0.165 


0.0444 


134. 


0.110 


0.0342 


31.6 


















(0.118) 


(0.0466) 


(54.3) 


O 7~» f~f O 


90 


400 


0. 


10 


0.582 


0.154 


53.0 


0.179 


0.134 


1553. 


















(0.178) 


(0.173) 


(25.0) 


O 7") 

d.P,S32 


140 


600 


0. 


10 


0.376 


0.0421 


114. 


0.098 


0.0183 


3.36 


















(0.102) 


(0.0339) 


(2.26) 




190 


onn 

800 


0. 


10 


0.176 


0.0187 


113. 


0.0636 


0.00574 


1.6 


















(0.0647) 


(0.0141) 


(1.00) 




1450 


onn 

300 


0. 


10 


0.000193 


0.0666 


150 


0.0433 


0.488 


1.84 


















(n r\A7n\ 

[U.U^i 1 U ) 


(^U.Ool ) 










n 

u. 


1 n 


0.00041 


0.0112 


134 


0.0490 


0.162 


0.611 


















(0.0490) 


(0.161) 


(0.606) 




2050 


600 


0. 


10 


0.000265 


0.00336 


125 


0.0207 


0.0585 


0.281 


















(0.0214) 


(0.0586) 


(0.280) 




200 


400 


0. 


30 


0.197 


0.0651 


120 


0.337 


1.10 


6.06 


















(0.350) 


(1.40) 


(4.17) 


^(/<?5)i 


220 


400 


0. 


40 


0.238 


0.0635 


121 


0.574 


19.1 


4.03 


















(0.620) 


(90.6) 


(3.97) 


AD(fg5) 2 


400 


900 


0. 


40 


0.0524 


0.0153 


113 


0.0516 


10.4 


0.814 


















(0.0528) 


(0.823) 


(0.828) 


SPS4 


400 


300 


0. 


50 


0.0424 


0.0685 


130 


5.69 


14.7 


7.5 


















(7.44) 


(9.06) 


(7.38) 


R0SZ2 


1000 


1000 


0. 


50 


0.014 


0.00118 


114 


1.41 


0.0341 


0.593 


















(0.640) 


(0.0295) 


(0.597) 
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Table 2: Non-universal mSUGRA Models: Input parameters are given at the grand 
scale using dimensions in GeV and the convention M 2 > 0. In all cases ji > 0. 





AD(fg9) 


5P56i 


SPS6 2 


SPS63 


CDG200 


CDG24 


CDG 75 


CDG011 


Mi 


420 


480 


720 


960 


800 


160 


-400 


424 


M 2 


420 


300 


450 


600 


160 


960 


240 


200 


M 3 


420 


300 


450 


600 


80 


-320 


80 


40 


m 


600 


150 


225 


300 


1400 


1400 


1400 


1400 


m H u 


600^ 


150 


225 


300 


1400 


1400 


1400 


1400 


A 


420 











1000 


1000 


1000 


1000 


tan (3 


40 


10 


10 


10 


50 


50 


50 


50 




Vijv(XiX" 


77) at Vij/c ~ 10 


i in units of 10 29 cm 3 


sec -1 






3.07 


0.482 


0.0973 


0.0767 


112 


0.000444 


0.00214 


72.2 




0.0102 


7.58 


1.36 


0.386 


7.56 


0.0339 


0.00488 


11.3 


X° 2 X°2 


14.0 


115 


133 


133 


1.64 


9.66 


106 


1.1 




Vij<7(XiX" 


99) at v i3 


/c ~ 10 - 


i in units of 10 29 cm 3 


sec -1 




xixi 


18.3 


1.6 


0.672 


0.376 


1356 


0.118 


2.80 


13.6 


x\xl 


1.34 


0.986 


0.514 


0.427 


19.8 


47.6 


2.82 


0.425 


xix°2 


1.19 


5.61 


12.7 


3529 


4.77 


7.72 


4.70 


0.0286 




% 


v(xW - 


> gg) at Vij/c ~ 0.5 


in units of 10 29 cm 3 sec 1 




xU 


20.0 


2.00 


0.811 


0.455 


565 


0.131 


2.84 


11.9 


x\xl 


1.39 


1.62 


1.42 


2.39 


19.0 


68.5 


2.50 


0.407 


X°2X°2 


1.67 


9.42 


2418 


5.17 


4.86 


7.79 


8.20 


0.0291 
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Table 3: GMSB and AMSB models. Input parameters are at the big scale, and 
v ij a (XiX°j ~^ 77? 99) predictions are at %/c ~ 1CT 3 (%/c ~ 0.5). 



Parameters (Dimensions in 


GeV) 




Vij(?(XiX 


?-77, 


99) in 10 


29 cm 3 sec 1 


n > 


77 


99 








GMSB 


models; 


n e = n q = 


■ N mes - 


= 3 






Model 




Msusy 


tan (5 


xixi 






Xi'Xi' 


xU» 2 




SPS7! 


80000 


40000 


15 


2.12 


0.225 


36.6 


2.01 


48.6 


7.53 
















(2.40) 


(20.9) 


(6.87) 


SPS7 2 


120000 


60000 


15 


1.15 


0.124 


25.5 


0.409 


7.34 


2.61 
















(0.456) 


(5.03) 


(2.43) 


SPS7 3 


160000 


80000 


15 


0.686 


0.0545 


17.4 


0.146 


2.87 


1.04 
















(0.180) 


(2.14) 


(0.991) 


SPS8! 


200000 


100000 


15 


0.416 


0.121 


131 


0.181 


0.994 


4.78 
















(0.192) 


(0.947) 


(8.30) 


SPS8 2 


400000 


200000 


15 


0.119 


0.00834 


121 


0.0299 


0.0457 


7.13 
















(0.0298) 


(0.0445) 


(26.9) 


AMSB models 


Model 


m 


m aux 


tan/3 


X1X1 


ft® 


X^ 


xW 


xix° 2 


X%X% 


SPS9! 


450 


60000 


10 


170 


0.0355 


0.487 


0.229 


0.0237 


0.315 
















(0.233) 


(0.0240) 


(0.336) 


SPS9 2 


675 


90000 


10 


148 


0.0149 


0.217 


0.110 


0.00364 


0.179 
















(0.112) 


(0.00374) 


(0.185) 


SPS9 3 


900 


120000 


10 


135 


0.00822 


0.122 


0.0652 


0.000984 


0.116 
















(0.0659) 


(0.00102) 


(0.117) 
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